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ABSTRACT 

A  computer-controlled  data  collection  system  was  designed  and 
implemented  to  analyze  finite  amplitude  acoustic  standing  waves  in  a 
rectangular  air-filled  cavity  with  linear  and  wedge-shaped  boundary 
perturbations.  The  response  of  a  nonlinear  standing  wave  in  per¬ 
turbed  and  unperturbed  cavities  was  studied.  The  experimental 
results  showed  that  the  computer-obtained  data  was  in  qualitative 
agreement  with  that  previously  obtained  by  manual  procedures.  Quan¬ 
titative  comparisons  with  theory  could  not  be  made  because  the  prop¬ 
erties  of  the  cavity  did  not  match  those  for  which  theoretical 
predictions  are  available. 
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A  one-dimensional  mathematical  model  for  finite-amplitude 
standing  waves  in  a  real  cavity  was  studied  by  Coppens  and  Sanders 
[Refs.  1.  2]  in  1968.  In  1974,  this  model  was  extended  to  three 
dimensions.  Losses  and  resonance  frequencies,  inadequately 
described  by  the  classical  theory  of  wall  losses,  are  obtained  for  the 
real  cavity  from  measurements  obtained  in  the  linear  acoustic  regime. 
The  research  done  in  1974  showed  that  when  a  cavity  is  driven  near 
the  resonance  of  a  standing  wave  the  nonlinear  distortion  primarily 
involves  only  the  normal  modes  of  the  cavity  that  are  “familially 
related”  to  the  driven  mode.  Lane  [Ref.  31.  DeVall  [Ref.  4],  Slocum 
[Ref.  51.  and  Kuntsal  [Ref.  6]  examined  the  problem  and  showed  that 
experimental  results  are  in  excellent  agreement  with  the  theory 
except  when  a  degeneracy  exists  in  the  cavity. 

The  purpose  of  the  research  reported  in  this  thesis  was  to  design 
and  implement  a  computer-assisted  data  collection  and  analysis  sys¬ 
tem  for  measuring  the  finite-amplitude  distortion  in  a  perturbed 
cavity  with  degeneracies. 


When  t±ie  sound  intensity  in  an  acoustic  process  becomes  so  high 
that  the  principle  of  superposition  no  longer  applies,  we  enter  the 
domain  of  nonlinear  acoustics.  The  limits  of  linear  approximations  in 
the  equation  of  force,  mass  continuity,  and  state  are  exceeded  and  the 
resultant  vibrational  disturbance  at  any  point  in  the  medium  will  no 
longer  be  equal  to  the  sum  of  its  individual  components.  The  essence 
of  the  nonlinear  acoustic  mechanism  is  that  an  intense  wave  perturbs 
the  medium  in  which  it  exists  and  this  perturbation  alters  the  natural 
order  governing  the  wave’s  behavior.  The  result  is  that  segments  of 
the  disturbance  in  a  warmer  compression  travel  faster  than  those  in  a 
cooler  rarefaction.  When  this  happens,  the  acoustical  wave  alters  its 
shape,  steepening  as  it  travels  along.  The  alteration  is  infinitesimal  for 
w'eak  disturbances  and  gradually  increases  in  importance  with 
increases  in  both  amplitude  and  frequency  of  sound.  [Ref.  2] 

In  1968,  Kirchoff  [Ref.  7]  studied  the  distortion  of  intense  acous¬ 
tic  waves  in  an  unbounded  fluid.  This  study  was  extended  by  Lamb 
[Ref.  8],  Fay  [Ref.  9],  and  by  Keller  [Ref.  10). 

Keck  and  Beyer  [Ref.  1 1  j  approached  the  problem  of  unbounded 
propagation  by  using  perturbation  methods  and  Fourier  series 
representations  of  the  waveform.  The  model  of  Coppens  and  Sanders 
[Ref.  12]  dealt  wrtth  finite-amplitude  standing  waves  in  rigid  walled 
cavities.  They  extended  the  Keck-Beyer  perturbation  approach  to 
include  wall  losses  predicted  by  the  Rayleigh-Kirchoff  theory  [Ref.  7], 

1  1 


and  showed  excellent  agreement  with  experiments  conducted  in  a 
rigid-walled  tube  at  low  but  finite  levels  of  excitation.  The 
experimental  work  by  Beech  (Ref.  13)  and  by  Ruff  [Ref.  141  showed 
that  at  high  excitation  levels  differences  developed  between  theory 
and  experiment.  Experimental  work  by  DeVall  [Ref.  4J  and  by  Kilmer 
[Ref.  15)  showed  that  if  degeneracies  existed  the  model  failed  to 
account  for  the  experimentally  observed  excitation  of  “nonfamiliar 
modes. 

The  effects  of  various  geometrical  boundary  perturbations  on  a 

rigid-walled  cavity  investigated  by  Aydin  [Ref.  16)  showed  that  the 

standing  waves  that  exist  in  an  ideal  cavity  must  be  corrected  when 

the  boundaries  are  irregular.  The  model  of  Coppens  and  Sanders  is 

based  on  a  three-dimensional,  nonlinear  wave  equation  with  a 

dissipative  term  describing  absorptive  losses  encountered  by  standing 

waves  in  rigid  walled  cavities.  As  stated  in  [Ref.  121, 

The  standing  waves  which  can  be  excited  in  any  real  cavity  deviate 
from  predictions  of  the  linear  wave  equation  with  ideal  boundary 
conditions  for  the  following  reasons: 

a  The  presence  of  boundary-layer  losses  at  the  cavity  surfaces 
yields  a  dispersive  contribution  to  the  wave  equation. 

b.  Geometrical  irregularities  alter  the  effective  dimensions  of  the 
cavity. 

Both  of  these  mechanisms  can  be  treated  separately  as  long  as  the 
shift  in  frequency  is  so  small  that  the  actual  resonances  are  close  to 
the  theoretical  values  resulting  from  the  classical  model.  These  are 
treated  by  assuming  the  dimensions  are  exact,  and  the  apparent  phase 
speeds  are  determined  on  that  basis. 
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The  resonance  frequency  for  each  standing  wave  is  defined  as 


(Ref.  6] 


o)n  =  Cn  [(nx  kx)2  +  (ny  ky)2  +  (nz  kz)2]1/2  ( 1 ) 

.  i  n  Jt  .  m  .  lx  , 
where  kx  =  ky  *  kz  =  and  n,  m.  and  1  are  integers.  Cn  is 

the  apparent  phase  speed  appropriate  for  that  frequency. 

The  force  equation,  in  Eulerian  coordinates,  can  be  written  in  the 
form 


f-+ta  -m  H-i-vp-ix* 


where  the  operator  X  describes  those  physical  processes  leading  to 
absorption  and  dispersion.  If  we  ignore  rotational  effects  so  that 
Tj  =  V0.  restrict  waves  to  rectangular  cavities,  and  neglect  terms  of 
orders  higher  than  M2.  M(a/k).  and  (a/k)2.  the  result  is  a 
quadratically  nonlinear  wave  equation  [Ref.  12]. 


The  standing  pressure  wave  in  a  rectangular  cavity  has  the  form 


[Ref.  12] 


p  =  coskx  x  coskyy  cosk2z  coscot 


...  >  ^  ***  W  ^  A  A  A  •*  A 


If  the  cavity  is  being  driven  near  resonance,  the  nonresonant  terms 
contribute  almost  nothing  with  respect  to  the  resonant  terms.  So  the 
standing  wave  pressure  has  the  form  of 

oo 

P  =  X  Pn  (5) 

n=  1 


where  each  Pn  is  of  frequency  nco  and  has  the  form  [Ref.  12]. 


Pn 


=  MRn  cosnkxx  cosnkyy  cosnkzz  sin(ncot  +  Qn) 


(6) 


These  pns  constitute  those  normal  modes  familially  related  to  pi. 

The  nonlinear,  coupled,  transcendental  equations  which  must  be 
solved  for  the  Rns  are 


cosl 

sinj 


*  (Qj  -  Qn-j) 


cosl 

sinj 


*  (Qn-j  -  Qj) 


for  all  n  >  1  (7) 


R  =  1  for  n=l 
where  Qn  =  Arctan(-Fn) 


Fn  =  Q  [(nco)2  -  (con)2]/(nco)2  =  2gn 


for 


nco  -  ton 
con 


«  1 


(8) 


Qn  - 


(On 


COnu  ~  G)nL 


(9) 


Total  pressure  near  resonance  for  perturbed  boundaries  is  [Ref.  6] 
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(10) 
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3 

*a 
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p  =  (n,  m.  1/g),0)  - 


oo 

vTw  X  aml  An  (~l)n  * 

*  SS i 


Snml  (n,m.l/co,0  +  Cnml) 


where  Snmi  -  first-order  perturbation  correction  and  onmi  angle 
which  is  defined  as  (Ref.  6]: 


-M 


A  =  Amp.  of  perturbation 


If  the  cavity  has  perturbed  walls,  the  solution  will  be  in  terms  of  a 
summation  of  the  classical  linear  solution  for  ideal  boundaries  plus 
perturbation  correction  terms  due  to  the  irregular  boundary: 

p  =  po  +  e  p’  +  e 2  p"  +  ...  (11) 

Since  the  magnitude  of  the  boundary  perturbation  is  kept  small,  sec¬ 
ond  and  higher  terms  in  e  can  be  considered  insignificant,  so  that 

p  =  po  +  €  p’  (12) 

and  p  must  satisfy  the  following  conditions. 

—  L  P  =  0  (13) 


yj 


j 

I 


53 


a 


$ 


Vp  it  =  0 


at  x  =  0,  Lx  [1  +  e  f(y.z)] 

y  =  0.  Ly 

z  =  0,  Lz 


where  it.  the  local  normal  to  the  real  surface,  is  obtained  by  taking  the 
gradient  of  the  equation  for  the  boundary  [Ref.  1), 

n=V{x-Lx(l+e  f(y,z)]}  ( 1 5) 

Thus,  to  the  first  order  in  e , 

A  T  3f(y.z)  A  .  af(y.z)  A  _ 

n  =  x  -  e  Lx  jjy  y  -  e  Lx  ^  z  (16) 

and  when  equation  (16)  is  used  in  equation  (15),  the  result  is 


'§£  af(y-z)  3p  _T  9f(y.z)  9pl 
LaX"6Lx  ay  dy~eLy  dz  dz\x  =  Lx[1  + 


|_dx  ay  ay  az  c/zj*  -  [i  +  ef[y>z)l 

Applying  a  Taylor-series  expansion  for  p  evaluated  at  the  real 
boundary  Lx  [1  +  ef(y,z))  produces 

pi  =  p  I  +  I  e  Lxf(y,z) 

x  =  Lx  (1  +  efly.zj]  x  =  Lx  °xx  =  Lx 


+  I  .  [eLxf(y,z)l2  +  ... 

OA  x  =  Lx 


Substituting  equation  (12)  into  the  RHS  of  equation  (18).  taking 
the  partial  derivative  with  respect  to  x  on  both  sides  and  keeping  the 
first-order  terms  in  e.  yields  (Ref.  1) 

H'  =  L,  [I  ♦  efly.zIl'H,1  =  4.  +  £  H'  •  L*  +  ^*U  €  W(y'Z’  +  ••  U9) 


Taking  the  partial  derivatives  with  respect  to  y  and  z  and  using 
exactly  the  same  procedure  gives  [Ref.  1] 

HU  II  ♦  eily.zl|=  HxU  +  *  HU  +  I^U  *  Wfly-Zl  +  -  1201 


H'=  Lx  II  *  eRy.z,f  Hx  Lx  +  S  HcL  Lx  +  ^xU  6  My'Z>  +  -  12  U 


Substituting  equations  (19),  (20).  and  (21)  into  equation  (17)  and 
keeping  the  first-order  terms  in  e  results  in 


sr_i  -  Lv  r_ f(v  z)  a2p.°  .  9f(y't)  9Po  +  8f(y-2)  aP°| 

^  ^  L  fly’z)  +  3y  3y+  3z  1TJX  =  ] 


The  RHS  of  equation  (22)  can  be  represented  as  a  Fourier  series  in 
cosines,  so  that  p'  can  be  expressed  as  a  summation  of  normal  modes. 
Hence, 


3pli  v  V  mit  Lk  .  „  .. 

1  =  U  U  ami  cos -f — y  cos  f— z  cos  (cot  +  9) 


-x  =  Lx  m=0  1=0 


(23) 


S 


•S3 


vA  ,v  v 


a 


1 

te 


i 

« 'l 

3 

Isb 


7&-I  _  =  £  £  ami  (O.m.l/w.0) 

UAv  =  1_A 


0Xx  =  Lx  ^oiTo 
where 

l  r7  r* 

aoo  -  t  t  J  J  (Gj  dy  d2 
0  0 


2  '?Y 


amO  =  -QQ  j  J  [G]cos^ydyd2 


2 


aoi  =  Or/  J  [G|  cos  -T—  dy  dj 


.  Lv  Lz 

4  fy  7  f  mrt  Ln  ,  m*0 

J  IGlcos-j^ycosj^zdyd^  n,0 


and  G  is  defined  as 


G#l 

ax*=4. 


In  order  to  find  the  contribution  to  p1  from  each  of  the  terms,  it  is 
stated  [Ref.  16]  that,  for  a  cavity  forced  by  a  dynamic  boundary  condi¬ 
tion  at  a  boundary, 

C^=lP=0  (30) 

and  the  dynamic  boundary  condition  from  the  (m.l)th  term  is 


■jrr- 1  =  A  cos  -r—  y  cos  7-  z  e«“t+Q> 

dx  x  =  Lx  h  Lz 


then 


*  1  v  a  /t  tm  m  n  mjc  Ltt 

A  ToTy —  ^ An  (_1)  Snmi  cos  l^x  cos  ~y cos  ~q 

p  Loc 
VCoJ 


.  gUcot+Q+OnmlJ 


where 


1  if  n  =  0 

An  "  [2  if  n  =  1.2.3.. 


Applying  this  solution  to  equation  (24)  gives  the  complete  solution 
for  the  first-order  perturbation  (Ref.  61 


p  =  (n.m,l/co.0)  - 


00 

w  \o  ^  aml  An  (—  1  )n  Snml 
v-o  J  m=o 


(n.m.l/w.0+Onml) 


As  it  is  stated  (Ref.  16],  if  this  is  near  a  resonance.  co-tonml.  then 
this  term  will  dominate  the  summation  and  all  other  non-degenerate 
terms  can  be  omitted.  Consequently,  equation  (33)  becomes 

p  =  (n.m.l/co.e)  — aml  An  (-l)n  Snmi  (n.m.l/G>,0+cnml)  •••  (34) 
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A.  COMPUTER-CONTROLLED  APPARATUS 

The  rectangular  cavity  used  in  this  research  was  designed  and 
constructed  in  1977  by  Kuntsal  (Ref.  61  (Figure  1).  One  wall  of  the 
cavity  could  be  moved  and  secured.  The  interior  of  the  cavity  was 
12.002  inches  long  and  2.502  inches  high,  while  the  width  could  be 
varied  from  approximately  5.50  inches  to  7.00  inches,  in  0.25-inch 
steps.  The  piston,  located  in  a  2.25-inch  diameter  port  in  the  floor  of 
the  cavity,  was  supported  by  an  O-ring  to  produce  a  seal  between  the 
piston  and  the  port.  Three  1 /4-inch  diameter  microphone  ports 
were  used.  A  condenser  microphone  was  mounted  in  a  0.24-inch 
aluminum  case  and  an  O-ring  and  silicon  grease  were  used  between 
this  case  and  the  ports  for  sealing. 

Modal  designation  for  the  lowest  three  modes  is  shown  in 
Figure  2.  The  (020)  and  (100)  modes  are  nearly  degenerate.  The 
(020)  mode  is  a  one-dimensional  standing  wave  with  nodes  3  and  9 
inches  along  the  long  wall  of  the  cavity  and  the  (100)  mode  is  a 
standing  wave  with  one  node  in  the  middle  of  the  short  wall.  When 
the  microphone  is  placed  at  position  A,  the  (100)  mode  is  sensed 
without  contamination  from  the  (020)  mode.  In  position  C,  the  (020) 
mode  is  sensed  without  contamination  from  the  (100)  mode.  A  block 
diagram  of  the  apparatus  is  shown  in  Figure  3. 

The  piston  was  driven  by  an  M-B  Electronics  model  EA1500 
exciter,  which  in  turn  was  driven  by  an  APS  Electronic  model  114 


(0.1.0) 


(0,2.0) 

(b) 


(1.0.0) 


Figure  2 

(a)  Cavity  Orientation 

(b)  Modal  Designation 


power  amplifier.  The  driving  signal  was  produced  by  an  HP  3325A 
Synthesizer/Function  generator. 

The  piston  movement  was  continuously  analyzed  by  an  Endevco 
model  2215  Accelerometer  mounted  within  the  piston.  The  output 
voltage  of  the  microphone  and  output  voltage  of  the  accelerometer 
were  measured  by  an  HP  3678  multimeter.  An  HP  3562A  Dynamic 
Signal  analyzer  was  used  to  measure  the  harmonic  distortion.  A  model 
5204  Lock-in  Analyzer  was  used  to  measure  the  resonance  frequen¬ 
cies  at  small  signal  amplitude. 

The  sound  pressure  in  the  cavity  was  sensed  by  a  Bruel  and  Kjaer 
4136  condenser  microphone  with  matching  preamplifier  2807.  The 
output  of  this  preamplifier  went  to  three  pieces  of  equipment:  (1)  an 
HP  3478  multimeter  to  measure  the  voltage  level.  (2)  an  HP  3562A 
Dynamic  Signal  Analyzer  to  measure  the  harmonics,  and  (3)  a  Cos 
5060A  oscilloscope  to  monitor  the  undesired  distortion  of  the  accel¬ 
erometer  signal. 

The  HP  3562A  Dynamic  Signal  Analyzer  and  HP  3325A  Synthe¬ 
sizer/Function  generator  were  controlled  by  an  HP  87  computer.  Two 
types  of  geometrical  perturbations  were  used: 

a  The  long  wall  was  machined  at  a  small  angle  (as  shown  in 
Figure  4). 

b.  The  long  wall  was  machined  as  shown  in  Figure  5. 

The  effect  of  these  geometrical  perturbations  on  the  resonance 
frequency  of  the  (100)  mode  could  be  estimated  by  calculating  an 
“effective"  width  Le  =  L  -  AL  of  the  perturbed  cavity  where  AL  =  V/A 
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Figure  4 


Linearly  Perturbed  Wall 

.25' 


Figure  5 

Wedged  Perturbation 

and  V  and  A  are  the  volume  and  cross-sectional  area  of  the 
perturbation.  [Ref.  1] 


B.  STRENGTH  PARAMETER  (STRM) 

Strength  of  the  finite  amplitude  interaction  is  characterized  by 
the  strength  parameter  STRM  [Ref.  6]. 

STRM  =  MPQi  where  Qi  is  the  quality  factor  of  the  driven  mode 


where  Sm  is  the  microphone  sensitivity. 


The  microphone  output  voltage  v*  is  read  on  an  HP  3562A 
Dynamic  Signal  Analyzer  set  at  the  fundamental  harmonic.  Micro¬ 
phone  sensitivity  was  obtained  with  a  Bruel  and  Kjaer  model  4220 
piston  phone. 

SM  =  U-56  ±  0.03)  x  10-3  volt/(N/M2) 

By  using  p  =  1.293  kg/m3  and  C  =  345  m/sec,  we  found 

STRM  =  7.07  x  10-3  x  Vx  g.  (35) 

This  is  the  same  sensitivity  as  is  measured  by  Kuntsal  [Ref.  6]. 

C.  FREQUENCY  PARAMETER  (FP) 

The  frequency  parameter  is  defined  by 


FP  =  2QL  (co -coi)/coi  (36) 

where  f  is  a  driving  resonance  frequency  and  fi  is  a  resonance  fre- 

Acon 

quency  of  (010)  mode.  It  is  assumed  that  — —  «  1.  FP  indicates 

u>n 

where  the  driving  frequency  is  with  respect  to  the  fundamental  reso¬ 
nance  frequency  of  the  cavity.  As  I  FP  I  increases,  the  nonlinear  effects 
will  in  general  become  smaller  for  a  given  strength  parameter. 


D.  HARMONICITY  COEFFICIENT  (E(n)) 

The  harmonicity  coefficient  is  defined  as  [Ref.  6) 


E(n)  = 


con  -  ncoi 
ncoi 


(37) 
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E(n)  characterizes  how  well  the  modes  of  a  given  family  are  tuned. 

E.  COMPUTER  PROGRAMS 

Two  computer  programs  were  developed.  The  first  program  finds 
the  resonance  frequencies  and  calculates  the  quality  factors  (Figure  6). 
The  second  program  drives  the  Dynamic  Signal  Analyzer  HP  3562A  by 
a  HP  3325A  Synthesizer/Function  generator  and  finds  and  calculates 
harmonic  content  of  the  microphone  output  at  each  driving  frequency 
(Figure  7). 
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Set  freq  *-  FR 


False  —  - 

Q55<FR<  1 700> - Stop 


Read  “FR1" 
Read  “FR2* 
Read  -FR3* 


CALCULATE  'O' 


Print  MAX(FR) 

Print  HaLfpower  (FR) 
Print  <Q) 


Figure  6 

’9”  Calculation  Flow  Diagram 
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Power  of  the  Harmonic  Content  Calculation  Flow  Diagram 


28 


-'.v.v 


'  "/V  V 


A.  SYSTEM  ALIGNMENT 
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The  primary  factor  in  setting  the  lower  limit  to  the  measurement 
of  harmonic  content  was  the  amount  of  distortion  injected  into  the 
cavity  by  the  piston.  Misalignment  of  the  piston  within  its  port  causes 
distortion  of  the  piston  motion  which  was  detected  by  analyzing  the 
harmonic  content  of  the  accelerometer  output  voltage.  (With  the  pis¬ 
ton  removed  from  its  port,  distortion  levels  were  of  the  same  order  of 
magnitude  as  the  general  electronic  noise  level  of  the  accelerometer 
system  when  at  rest.)  Before  each  run.  the  piston  was  aligned  within 
its  port  so  that  the  accelerometer  output  voltage  contained  less  than 
0.1  percent  second  harmonic  and  0.2  percent  third  harmonic. 

After  testing  various  configurations,  the  one  selected  for  ease  of 
alignment  consisted  of  the  cavity  mounted  on  a  tripod  support  directly 
above  the  exciter,  which  was  supported  by  its  own  base  in  an  upright 
position.  Both  the  cavity  tripod  and  the  exciter  rested  on  a  5/ 16-inch 
pad  of  stiff  acoustic  rubber.  Once  the  system  was  aligned,  only  an 
occasional  minor  adjustment  was  necessary  to  maintain  low  piston 
distortion. 

B.  AUTOMATIC  DATA  COLLECTION 

Since  the  resonance  frequencies  of  the  cavity  were  observed  to 
drift  with  time,  the  system  was  allowed  to  warm  up  for  at  least  three 
hours  prior  to  data  collection.  The  harmonic  content  of  the 
accelerometer  output  was  analyzed  and  the  piston  was  adjusted  until 


the  second  harmonic  of  the  accelerometer  output  was  at  least  50  dB 
below  that  of  the  fundamental  frequency.  Figure  8  shows  the  drift  in  fr 
is  approximately  linear  in  time. 

Data  was  collected  using  two  measurement  procedures  in  the 
computer  programs:  (1)  a  pre-  and  post-run  infinitesimal  amplitude 
measurement:  and  (2)  the  finite  amplitude  measurement. 

1.  Pre-Run  and  Post-Run  Infinitesimal  Amplitude 
Measurements 

During  the  pre-run  and  post-run  infinitesimal  amplitude 
measurement,  the  piston  was  driven  80  mV  (rms).  Because  the  (100) 
and  (020)  modes  were  nearly  degenerate,  measurements  of  the  prop¬ 
erties  of  the  (100)  mode  were  made  with  the  microphone  in  port  C 
while  measurements  of  the  (020)  mode  were  made  with  the  micro¬ 
phone  in  port  B.  The  resonant  frequency  of  the  fundamental  and  the 
Q  factor  first  few  (O.n.O)  in  the  modes  were  determined  from 

fr  =  i'}L^k  (38) 


9  = 


fr 

fu-fL 


(39) 


and  the  harmonicity  coefficients  were  then  calculated  from 

p  fn~nfi 
En  ”  nfi 


(40) 
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Time  (ks) 


Figure  8 

Drift  in  Frequency 


where  fu  and  fi,  are  the  frequencies  of  the  upper  and  lower  half-power 
points,  and  fn  is  the  resonance  frequency  of  the  nth  member  of  the 
(010)  family.  Figure  9  shows  the  precision  of  frequency 


measurements. 

2.  The  Finite-Amplitude  Measurement 

All  finite-amplitude  measurements  were  made  at  port  B. 
Throughout  this  measurement.  Vj  had  to  be  kept  constant  by  the 
computer  program  (Appendix  C).  which  required  adjustment  of  the 
driving  voltage  applied  to  the  piston.  The  frequency  was  increased  by 
the  computer  in  0.5-Hz  steps  through  ±5  Hz  about  the  resonance 
frequency  of  a  fundamental  mode.  At  each  driving  frequency,  the 
harmonic  content  of  the  microphone  output,  usually  up  to  the  fourth, 
was  measured  with  the  HP  3562A  Dynamic  Signal  analyzer. 

To  determine  the  exact  value  of  the  frequency  parameter,  the 
value  of  fr  at  the  instant  of  the  measurement  of  Vn  was  determined 
from  the  first  computer  program. 
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V.  RESULTS 


In  this  section,  the  theoretical  predictions  and  experimental 
results  obtained  by  Kuntsal  (Ref.  6]  will  be  qualitatively  compared  to 
the  experimental  results  obtained  by  using  the  computer-controlled 
equipment. 

The  information  on  the  empirical  losses  and  resonance  frequen¬ 
cies  is  contained  in  the  Q's  and  E’s.  It  is  important  to  note  that, 
because  V]  is  always  kept  constant  as  a  function  of  frequency,  the  sec¬ 
ond  harmonic  distortion  is  maximized  for  co  =  a>2/2.  The  point  where 
this  occurs  for  each  V^/Vi  curve  is  indicated  by  an  arrow  with  the 
label  Fq20-  At  that  point,  the  value  of  the  frequency  parameter  is 

F020  =  2  Qoio  E020 

Similarly. 

F030  =  2  Qoio  E030 

is  indicated  on  each  graph.  The  arrow  labelled  Fioo  indicates  the 
position  where  the  nearly  degenerate  (1,0,0)  mode  is  resonant,  and 
the  value  of  Fioo  is 

F ioo  =  2  Qoio  ElOO 

The  theoretical  predictions  for  unperturbed  cavities  made  by 
Kuntsal  (Ref.  6]  were  compared  to  the  experimental  results  as  seen  in 


Figure  A.  It  is  clearly  seen  that  the  theory  and  experiment  are  qualita¬ 
tively  in  agreement. 

At  Figure  9.  fi  and  f2  are  plotted  as  a  function  of  time.  This  figure 
showed  two  things:  (1)  changes  of  resonance  frequencies  of  the  family 
modes  are  proportional  to  the  frequency,  as  would  be  expected  if  the 
long  dimension  of  the  cavity  changes  with  time,  most  likely  due  to 
temperature  changes:  and  (2)  the  fluctuations  of  resonance  frequen¬ 
cies  from  measurement  to  measurement  indicate  a  precision  of  ±  0.02 
Hz  in  determining  the  resonance  frequencies  of  the  lowrest  mode.  By 
using  the  value  of  ±  0.02  Hz,  expected  uncertainties  in  the  E's  are 

fn  ±  0.02  -  n(fi  ±  0.02) 

^Eexpected  = 

Calculated  AEexpected  and  measured  AE's  between  pre-  and  post-run 
are  tabulated  in  Table  4.  The  difference  of  the  pre-run  and  post-run 
E's  for  the  family  modes  in  the  unperturbed  cavity  are  within  the 
expected  precision.  For  the  100  mode  (non-family  mode),  this  differ¬ 
ence  is  much  greater  than  expected,  indicating  that  the  cavity  may  be 
changing  its  dimensions  unequally  in  all  directions. 

For  a  linear  perturbation  (Figure  B)  and  for  the  wedge  perturba¬ 
tion  (Figure  C).  the  harmonic  content  in  Reference  6  is  in  qualitative 
agreement  with  experiment,  considering  that  lower  Q's  and  different 
E's  were  observed  compared  to  Kuntsal's  measurement  [Ref.  6]. 
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The  experimental  apparatus,  procedures,  and  computer  programs 
are  capable  of  providing  data  sufficiently  precise  to  test  the  prediction 
of  the  theory  in  the  absence  and  presence  of  geometrical  perturbation. 

Because  the  strength  parameters  (STRM)  and  quality  factors  (Q) 
calculated  in  .the  experiment  did  not  match  the  STRM  and  Q  in  Refer¬ 
ence  6,  quantitative  agreement  could  not  be  tested.  However,  the 
experimental  results  show  that  a  geometrical  perturbation  alters  the 
finite  amplitude  behavior  of  the  cavity,  and  that  the  nature  of  these 
changes  is  in  qualitative  agreement  with  the  predictions  of  the  theory. 

The  quantitative  disagreements  between  the  theory  and  experi¬ 
ment  for  the  perturbed  cavities  may  be  because  the  resonance  fre¬ 
quency  of  the  non-family  mode  is  not  tracking  the  family  modes. 

In  order  to  avoid  disturbing  the  position  of  the  piston  in  the  port 
of  the  cavity,  the  microphone  has  to  be  shifted  v^ry  carefully  from  one 
port  to  another  between  the  "pauses'  as  the  program  goes  on. 


Average  Values  for  9s  and  Es 


Mode:  010 
STRM:  .370 


9,  = 
©2  = 
©3  “ 

©ioo  ~ 


187.52 

288.13 

350.42 

291.40 


E,  =  0.0 

E2  =  -49.02  x  10-4 

E3  =  -23.25  x  10-4 

E100  =  -51.34  x  104 


APPENDIX  B 

tables 

TABLE  1 

UNPERTURBED  CAVITY 

(SEE  FIGURE  A) 


Pre-Run 


Mode 

fu 

fL 

fr 

9 

E 

010 

565.801 

562.875 

564.351 

195.40 

0.0 

020 

1131.840 

1127.791 

1 129.812 

278.69 

98.16xl0-5 

030 

1698.271 

1693.524 

1695.901 

356.51 

16.82xl0-4 

100 

1131.872 

1127.590 

1 129.731 

264.20 

91.17x10-5 

Finite-Amplitude  Run 


Vi  =  -11  dB  Mode:  010  STRM  =  .389 


f  (Hz) 

V2  (dB) 

V3  (dB) 

V4  (dB) 

V2/Vi 

V3/V1 

V4/V  ! 

FP 

558.5 

-50.0 

-70.5 

-70.5 

0.0113 

0.0014 

0.0012 

-4.19 

559.0 

-47.9 

-72.0 

-69.0 

0.0143 

0.0013 

0.0011 

-3.85 

559.5 

-46.9 

-73.8 

-67.8 

0.0160 

- 

0.0011 

-3.50 

560.0 

-46.0 

-72.6 

-67.5 

0.0179 

- 

0.0010 

-3.15 

560.5 

-44.8 

-71.2 

-69.7 

0.0204 

0.0011 

0.0010 

-2.81 

561.0 

-43.7 

-69.7 

-73.5 

0.0233 

0.0012 

— 

-2.47 

561.5 

-42.4 

-67.8 

-76.5 

0.0269 

0.0014 

— 

-2.12 

562.0 

-41.0 

-65.7 

-75.9 

0.0316 

0.0018 

- 

-1.78 

562.5 

-39.4 

-63.2 

-82.5 

0.0380 

0.0025 

— 

-1.43 

563.0 

-37.8 

-60.3 

-85.0 

0.0457 

0.0034 

- 

-1.09 

563.5 

-35.9 

-56.6 

-74.9 

0.0569 

0.0052 

- 

-0.74 

564.0 

-34.0 

-52.2 

-67.5 

0.0708 

0.0070 

0.0015 

-0.40 

564.5 

-32.2 

-47.6 

-60.6 

0.0876 

0.0148 

0.0033 

-0.05 

565.0 

-31.1 

-45.0 

-54.5 

0.0994 

0.0200 

0.0067 

0.29 

565.5 

-31.6 

-43.9 

-52.9 

0.0933 

0.0226 

0.0070 

0.64 

566.0 

-33.3 

-46.7 

-56.0 

0.0772 

0.0164 

0.0056 

0.98 

566.5 

-35.3 

-51.0 

-63.2 

0.0610 

0.0100 

0.0025 

1.33 

567.0 

-37.6 

-55.5 

-70.8 

0.0484 

0.0060 

0.0010 

1.67 

567.5 

-39.8 

-60.0 

-73.6 

0.0365 

0.0036 

- 

2.02 

568.0 

-42.1 

-64.3 

-74.7 

0.0279 

0.0022 

- 

2.36 

568.5 

-44.1 

-68.5 

-74.2 

0.0221 

0.0013 

- 

2.71 

569.0 

-46.2 

-69.5 

- 

0.0174 

0.0011 

- 

3.05 
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TABLE  2 
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LINEAR  PERTURBATION  CONFIGURATION 

(SEE  FIGURE  B) 

Pre-Run 


Mode 

fu 

fL 

fr 

9 

E 

■ 

569.465 

1133.883 

1703.014 

1 131.650 

566.545 

1129.961 

1698.008 

1131.341 

565.653 

1132.306 

1700.594 

1 134.651 

193.72 

288.71 

339.72 
263.321 

0.0 

17.67xl0-4 

21.42xl0-4 

29.57x10-4 

Finite-Amplitude  Run 

Vx  =  -12  dB  Mode:  010  STRM  =  .355 


f  (Hz) 

V2  (dB) 

V3  (dB) 

V4  (dB) 

V2/Vi 

V3/V! 

V4/V ! 

FP 

561.0 

-43.0 

-68.7 

0.0251 

0.0013 

-3.19 

561.5 

-42.7 

-68.0 

- 

0.0260 

0.0014 

- 

-2.84 

562.0 

-42.3 

-66.9 

- 

0.0271 

0.0016 

- 

-2.50 

562.5 

-41.6 

-64.5 

- 

0.0296 

0.0021 

- 

-2.16 

563.0 

-40.3 

-62.7 

- 

0.0342 

0.0026 

- 

-1.82 

563.5 

-38.7 

-60.1 

- 

0.0411 

0.0035 

- 

-1.47 

564.0 

-37.3 

-56.2 

- 

0.0485 

0.0055 

- 

-1.13 

564.5 

-34.2 

-53.5 

- 

0.0691 

0.0075 

- 

-0.79 

565.0 

-33.8 

-51.0 

-65.0 

0.0721 

0.0100 

0.0020 

-0.44 

565.5 

-31.8 

-46.4 

-60.1 

0.0911 

0.0170 

0.0035 

-0.10 

566.0 

-30.9 

-43.0 

-54.6 

0.1012 

0.0251 

0.0066 

0.24 

566.5 

-29.6 

-40.3 

-50.0 

0.1161 

0.0344 

0.01 12 

0.58 

567.0 

-30.9 

-45.8 

-54.1 

0.1000 

0.0180 

0.0070 

0.92 

567.5 

-32.8 

-49.8 

-59.2 

0.0810 

0.0115 

0.0039 

1.26 

568.0 

-34.1 

-51.5 

-65.0 

0.0700 

0.0094 

0.0020 

1.61 

568.5 

-36.7 

-56.2 

-68.1 

0.0514 

0.0055 

0.0014 

1.95 

569.0 

-38.0 

-58.7 

- 

0.0441 

0.0041 

- 

2.29 

6' 


Post-Run 


Mode 

fu 

fL 

fr 

Q 

E 

M 

570.410 

567.456 

566.861 

191.90 

0.0 

1 135.067 

1131.116 

1133.984 

287.01 

23.10x10-4 

mm 

1705.321 

1700.340 

1701.110 

341.50 

39.98X10-4 

100 

1 136.740 

1132.416 

1135.756 

262.66 

17.94x10-4 
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TABLE  3 


WEDGE  PERTURBATION  CONFIGURATION 

(SEE  FIGURE  C) 


Pre-Run 


Mode 

fu 

fL 

fr 

9 

E 

B 

566.912 

1127.191 

1695.310 

1126.761 

563.864 

1123.281 

1690.471 

1122.911 

565.381 

1125.230 

1692.891 

1124.829 

185.80 

288.08 

349.27 

291.41 

0.0 

-48.92xl0-4 

-19.17x10-4 

-52.47x10-4 

Finite-Amplitude  Run 

Vi  =  -11  dB  Mode:  010  STRM  =  .370 


f  (Hz) 

V2  (dB) 

V3  (dB) 

V4  (dB) 

V2/V! 

v3/v  1 

V4/V! 

FP 

559.4 

-39.5 

-66.9 

0.0376 

0.0016 

-3.99 

559.9 

-37.9 

-65.3 

- 

0.0450 

0.0019 

- 

-3.66 

560.4 

-36.4 

-63.9 

_ 

0.0540 

0.0023 

- 

-3.33 

560.9 

-34.5 

-61.8 

- 

0.0672 

0.0029 

- 

-3.00 

561.4 

-32.4 

-59.1 

-78.0 

0.0856 

0.0039 

- 

-2.66 

561.9 

-30.3 

-55.7 

-72.7 

0.1090 

0.0058 

0.0008 

-2.33 

562.4 

-28.7 

-55.6 

-68.2 

0.1310 

0.0059 

0.0014 

-2.00 

562.9 

-28.9 

-50.8 

-65.3 

0.1271 

0.0103 

0.0017 

-1.67 

563.4 

-30.9 

-50.0 

-62.5 

0.1010 

0.0112 

0.0027 

-1.34 

563.9 

-32.3 

-49.7 

-60.3 

0.0810 

0.0116 

0.0034 

-1.01 

564.4 

-34.7 

-50.9 

-61.4 

0.0650 

0.0101 

0.0030 

-0.67 

564.9 

-36.3 

-54.3 

-67.0 

0.0550 

0.0068 

0.0016 

-0.34 

565.4 

-37.5 

-59.3 

-75.2 

0.0480 

0.0038 

- 

-0.01 

565.9 

-38.4 

-65.0 

- 

0.0430 

0.0020 

- 

0.32 

566.4 

-39.3 

-69.7 

- 

0.0391 

0.0012 

- 

0.65 

566.9 

-39.7 

-72.2 

- 

0.0370 

0.0009 

- 

0.98 

567.4 

-40.6 

-74.5 

- 

0.0280 

- 

- 

1.32 

567.9 

-40.7 

-75.4 

- 

0.0230 

- 

- 

1.65 

568.4 

-41.6 

-73.5 

- 

0.0210 

- 

- 

1.98 

568.9 

-44.0 

-71.5 

- 

0.0193 

- 

- 

2.31 

569.4 

-47.3 

-70.9 

- 

0.0151 

- 

- 

2.64 

569.9 

-49.1 

-71.1 

- 

0.0120 

- 

2.97 

12‘ 


Post-Run 


[  Mode 

- ?u - 

*L 

- E - 

9 

E 

H 

566.950 

1127.311 

1695.544 

1126.981 

563.963 

1123.411 

1690.731 

1123.126 

565.464 

1 125.371 
1693.144 
1125.051 

189.24 

288.18 

351.57 

291.39 

0.0 

-49. 12xl0-4 
-27.33xl0~4 
-50.2  lx  10-4 

TABLE  4 

E’S  DIFFERENCES  BETWEEN  RUNS 


Wedge  Perturbed 
Configuration 

AE2 

. 416xl0-4 

5.43x10-4 

0.2xl0-4 

AE3 

.08x10-4 

8.56x10-4 

8.16x10-4 

AE100 

1.7x10-4 

1 1.63x10-4 

2.26x10-4 

-^Eexpeeted 


.53x10-4 


APPENDIX  C 


COMPUTER  PROGRAMS 


This  program  drives  the  HP  3325A  Synthesizer/Function  genera¬ 
tor  between  562  Hz  and  1710  Hz.  finds  the  resonance  frequency,  and 
calculates  the  quality  factor  (Q). 


1 0  !  SUBROUTINE  Q 

14  SETTIME  0,87200 

1 5  FOR  INC-1  TO  30 
20  K— 1 

30  FOR  FREQ-562  TO  1 71 0  STEP  2 

50  IF  FREQ>550  AND  FREQ<570  THEN  GOTO  210 

70  IF  FREQ>570  THEN  GOTO  90 

90  FOR  FREQ-1 116  TO  1140  STEP  2 

110  IF  FREQ>1 1 15  AND  FREQ<1 140  THEN  GOTO  210 

130  IF  FRQ>1 140  THEN  GOTO  1791 

140  K-K+1 

150  FOR  FREQ-1680  TO  1710  STEP  2 

170  IF  FREQ>1670  AND  FREQd  710  THEN  GOTO  210 

190  IF  FREO>1710  THEN  GOTO  1800 

210  FI -FREQ 

230  F2-F1 +2 

250  F3-F1+4 

270  REMOTE  717,723,702 

290  OUTPUT  71 7  ;TR";F1  ;"HZ" 

310  WAIT  5000 

330  ENTER  723  ;AMP1 

350  OUTPUT  71 7  ;"FR";F2;"HZ 

370  WAIT  5000 

390  ENTER  723  ;  AMP2 

4 1 0  OUTPUT  71 7  ;"FR";F3;"HZ" 

430  WAIT  5000 

450  ENTER  723  ;  AMP3 

470  DISP  FI  .F2.F3.AMP1  ,AMP2,AMP3,K 

510  IF  AMP2<AMP1  OR  AMP2<AMP3  THEN  GOTO  590 

570  IF  AMP2>AMP1  OR  AMP2>AMP3  THEN  GOSUB  BINOM 

590  NEXT  FREQ 

610  !  SUBROUTINE  BINOM 


630  BINOM: 

650  P=F2 
670  S=2 

690  S=S/2 

710  OUTPUT  717  ;“FR";P;"HZ" 

720  WAIT  5000 

730  ENTER  723  ;  AP 

750  IF  P>1 71 0  THEN  GOTO  1 800 

770  PU=P+S 

790  OUTPUT  71 7  ;"FR";PU;"HZ" 

800  WAIT  5000 
810  ENTER  723  ;  AU 
830  PL=P-S 

850  OUTPUT  71 7  ;"FR";PL;"HZ" 

860  WAIT  5000 
870  ENTER  723  ;  AL 
890  DISP  P,PU,PL,AP,AU,AL 
910  IF  AU>AP  THEN  P=PU 
930  IF  AL>AP  THEN  P=PL 
950  IF  S>.008  THEN  GOTO  690 
970  GOSUB  HPWR 
990  RETURN 

1010  !  SUBROUTINE  HALFPOWER 
1030  HPWR: 

1050  DISPP.AP 
1060  DISP  TIME 
1065  PRINT  TIME 
1070  PRINT  "MAX- 
1090  PRINT  P,AP 

1100  IF  INC>1  AND  INC<30  THEN  GOTO  1 771  ELSE  GOTO  1110 
1110  PRINT  -REAL  HALFPOWER" 

1130  APH1=AP/2A.5 
1150  PRINT  APH1 
1 1 70  FOR  F-1  TO  2 
1190  PH=P+((-1  )AF) 

1210  OUTPUT  71 7  ;"FR";PH;"HZ" 

1230  WAIT  5000 
1250  ENTER  723  ; AH 
1270  DISP  PH, AH 

1290  IF  AH<APH1  THEN  GOTO  1350  ELSE  GOTO  1310 
1310  PH=PH+((-1  )AF) 

1330  GOTO  1210 
1350  S=.5 

1370  PHI=PH-(S*((-1  )AF)) 

1390  FOR  1=1  TO  10 

1410  OUTPUT  717  ;"FR";PHI;"HZ" 

1430  WAIT  5000 


1450  ENTER  723  ;  AHI 
1470  S=S/2 

1 490  IF  AHI>APH1  THEN  GOTO  1 570 
1510  FHI=PHI-(S*((-1  )AF)) 

1530  NEXT  I 

1550  GOTO  1610 

1570  PHI=PHI+(S*((-1  )AF)) 

1590  NEXT  I 

1610  PRINT  "HALFPOWER" 

1630  DISP  PHI, AHI 
1650  PRINT  PHI, AHI 
1670  IF  F=1  THEN  PHIR=PHI 
1960  NEXT  F 
1710  Q=P/(PHI-PHIR) 

1730  PRINT  "Q  VALUE- 
1750  PRINT  Q 

1770  DISP  Q 

1771  K=K+1 

1772  IF  K=0  THEN  GOTO  90 

1  773  !  AFTER  FIRST  PAUSE  SHIFT  THE  MICROPHONE  TO  THE  PORT  C  TO 
MEASURE  THE  FREQUENCY  OF  THE  DEGENERATE  MODE. 

1 775  !  AFTER  SECOND  PAUSE  SHIFT  THE  MICROPHONE  TO  THE  PORT  B 

1791  IF  K=1  THEN  PAUSE 

1792  IF  K=2  THEN  PAUSE 

1 793  IF  K=1  THEN  GOTO  90 

1794  IF  K=2  THEN  GOTO  140 
1797  RETURN 

1800  NEXT  INC 
1810  END 


This  program  drives  the  HP  3562A  Dynamic  SignalAnalyzer  by  HP 
3325A  Synthesizer/Function  generator  and  finds  and  computes  the 
harmonic  content  of  the  microphone  output  at  each  driving  frequency 
between  563  Hz  and  570  Hz. 


10  REMOTE  720,717 

1  5  SETTIME  0,87033 

20  OUTPUT  720  ;HPRST" 

30  OUTPUT  720  ;"B" 

40  OUTPUT  720  ;“FRS3KHZH 

50  OUTPUT  720  ;"PSPCM 

60  P=800 

70  OUTPUT  71 7  ;"AMP";P;"MV" 

80  OUTPUT  720  ;"LNRS" 

1 00  OUTPUT  720  ;"STBL" 


1 1 0  OUTPUT  720  ;"PSP2" 

1 20  OUTPUT  720  ;"KEY31  " 

1 25  !  DRIVING  THE  DYNAMIC  SIGNAL  ANALYZER  AND  READING  THE 

HARMONIC  CONTENT  OF  THE  DRIVING  FREQUENCY. 

1 30  FREQ-563 

1 40  OUTPUT  71 7  ;"FR";FREQ;"HZ" 

1 50  OUTPUT  720  SF563HZ  " 

160  WAIT  3000 
1  70  OUTPUT  720  ;"STRr 
1  80  WAIT  8000 
1 90  OUTPUT  720  ;"MKPK" 

200  OUTPUT  720  ;"RDMK" 

210  ENTER  720  ;  XI  ,YB 
220  DIL=YB 
230  DISPX1.DIL 

234  PRINT  TIME 

235  PRINT  XI, DIL 

240  FOR  FREQ=563  TO  570  STEP  .5 

250  OUTPUT  720  ;"FRS3KHZ" 

251  OUTPUT  720  ;"SF";FREQ;"HZ" 

255  OUTPUT  71 7  ;"FR";FREQ;"HZ" 

260  FOR  1=1  TO  4 

270  OUTPUT  71 7  ;"AMP";P;"MV" 

280  OUTPUT  720  ;"STRr 

290  WAIT  8000 

300  OUTPUT  720  ;"MKPK" 


50 


r.  r.  / 


% 

310 

OUTPUT  720  ;MRDMKM 

« 

320 

ENTER  720  ;  X,YB 

m 

324 

DISP  X,YB 

325 

PRINT  X,YB 

f 

326 

PRINT  TIME 

332 

IF  1=1  THEN  GOTO  333  ELSE  GOTO  350 

t 

333 

IF  FREQ>563  THEN  GOTO  335  ELSE  GOTO  350 

335 

GOSUB  BINOM 

350 

OUTPUT  720  ;"SF";X+50;"HZ" 

351 

OUTPUT  720  FRS.7KHZ" 

8 

355 

NEXT  1 

Jv 

356 

NEXT  FREQ 

& 

360 

!  SUBROUTINE  BINOM 

365 

!  THIS  SUBROUTINE  KEEPS  THE  FUNDAMENTAL  CONSTANT 

370 

BINOM: 

380 

RAT=1 0A(DILV1 0)/1 0A(YB/1 0) 

p.- 

fcj 

]OC 

385 

DISP  RAT.DIL.YB 

390 

IF  RAT>.999  AND  RAT<1.001  THEN  RETURN 

$ 

400 

P=P*RATA.2 

410 

OUTPUT  717  ;"AMP";P;"MV" 

1 

420 

OUTPUT  720  rsTRr 

430 

WAIT  8000 

i  v{i 

440 

OUTPUT  720  ;"MKPK" 

| 

450 

OUTPUT  720  ;"RDMK" 

460 

ENTER  720  ;  X.YB 

9 

470 

DISP  X.YB 

■j 

475 

PRINT  TIME 

a 

480 

PRINT  X.YB 

s 

490 

IF  FREQ>570  GOTO  660 

tv-: 

510 

GOTO  380 

* 

660 

END 
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